Resonance Damping in Ferromagnets and Ferroelectrics 



A. Widom 

Physics Department, Northeastern University, Boston, MA USA 

S. Sivasubramanian 
NSF Nanoscale Science & Engineering Center for High-rate Nanomanufacturing, 
Northeastern University, Boston MA USA 

C. Vittoria and S. Yoon 
Department of Electrical and Computer Engineering, Northeastern University, Boston, MA USA 

Y.N. Srivastava 

Physics Department and & INFN, University of Perugia, Perugia IT 

The phenomenological equations of motion for the relaxation of ordered phases of magnetized 
and polarized crystal phases can be developed in close analogy with one another. For the case of 
magnetized systems, the driving magnetic field intensity toward relaxation was developed by Gilbert. 
For the case of polarized systems, the driving electric field intensity toward relaxation was developed 
by Khalatnikov. The transport times for relaxation into thermal equilibrium can be attributed to 
viscous sound wave damping via magnetostriction for the magnetic case and electrostriction for the 
polarization case. 
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I. INTRODUCTION 

It has long been of interest to understand the close 
analogies between ordered electric polarized systems, e.g. 
ferroelectricity, and ordered magnetic systems, e.g. fer- 
romagnetism. At the microscopic level, the source of such 
ordering must depend on the nature of the electronic en- 
ergy spectra. The relaxation mechanism into thermal 
equilibrium state must be described by local electric field 
fluctuations for the electric polarization case and by mag- 
netic intensity fluctuations for the magnetization case; 
Speciflcally, the fleld fluctuations for each case 
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determine the relaxation time tensor for both cases via 
the fluctuation-dissipation formulai^'^— 



lim 



V JV 



dt. (2) 



We have unified the theories of relaxation in ordered po- 
larized systems and ordered magnetized systems via the 
Kubo transport time tensor in Eqs.([T]) and ©. 

The transport describing the relaxation of or- 
dered magnetization is the Landau-Lifshitz-Gilbert 
equation^"—. This equation has been of considerable 
recent interesliSi^— in describing ordered magnetic reso- 
nance phenomena^ 1^14. The equation describing the elec- 



tric relaxation of an ordered polarization is the Landau- 
Khalatinikov-Tani equationi^^— . This equation can be 
simply modeled^*"^^ with effective electrical circuits^^r^^. 
Information memory applications^^^^^ of such polarized 
system are of considerable recent intereslj^"— . 

The unification of the magnetic Gilbert-Landau- 
Lifshitz equations and the electric Landau-Khalatnikov- 
Tani equations via the relaxation time tensor depends 
on the notion of a nonequilibrium driving field. For the 
magnetic case, the driving magnetic intensity deter- 
mines the relaxation of the magnetization via the torque 
equation 



M = 7M X Hd 



(3) 



wherein 7 is the gyromagnetic ratio. For the electric case, 
the driving electric fleld determines the relaxation of 
the polarization via the equation of motion for an ion of 
charge ze 



mr = zeEj 



(4) 



The uniflcation of both forms of relaxation lies in the 
close analogy between the magnetic driving intensity 
and the electric driving field E^. 

In Sec[lT] the thermodynamics of ordered magnetized 
and polarized systems is reviewed. The notions of mag- 
netostriction and electrostriction are given a precise ther- 
modynamic definition. In Sec lIIIl the phenomenology of 
the relaxation equations are presented. The magnetic 
driving intensity and the electric driving field Ed are 
defined in terms of the relaxation time tensor Eq.Q. In 
Sec II VI we introduce the crystal viscosity tensor. From a 
Kubo formula viewpoint, the stress fiuctuation correlax- 
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ation 

1 

T^jki{r,r',t) = - {Aaki{r\-iX)Aa,j{r,t))dX, (5) 



determines the crystal viscosity 



Vijkl 



lim 



V Jv 



dt. 



(6) 



When the system is out of thermal equilibrium, the driv- 
ing magnetic intensity is 

wherein t are the relaxation time tensor transport co- 
efficients which determine the relaxation of the ordered 
magnetic system into a state of thermal equilibrium. 



For models of magnetic relaxation wherein acoustic heat- 
ing dominates via magnetostriction'^'^ and for models of 
electric relaxation wherein acoustic heating dominates 
via electrostriction, the relaxation time tensor in Eq.([2]) 
can be related to the viscosity tensor Eq.®. An inde- 
pendent microscopic derivation of viscosity induced re- 
laxation is given in Appendix!^ In the concluding SecfVl 
the sound wave absorption physics of the viscous damp- 
ing mechanism will be noted. 



II. THERMODYNAMICS 

Our purpose is to review the thermodynamic proper- 
ties of both magnetically ordered crystals and polariza- 
tion ordered crystals. The former is characterized by 
a remnant magnetization M for vanishing applied mag- 
netic intensity H ^ while the latter is characterized by 
a remnant polarization P for vanishing applied electric 
field E -> 0. 



A. Magnetically Ordered Crystals 

Let w be the enthalpy per unit volume. The funda- 
mental thermodynamic law determining the equations of 
state for magnetically ordered crystals is given by 



dw = Tds + H ■ dM -e:d(T. 



(7) 



wherein s is the entropy per unit volume, T is the temper- 
ature, e is the crystal strain and cr is the crystal stress. 
The magnetic adiabatic susceptibility is defined by 



If 



M 



(8) 



(9) 



denotes a unit vector in the direction of the magnetiza- 
tion, then the tensor Aijki describing adiabatic magne- 
tostriction coefficients may be defined as^i 



y^'^^iJ s.m 



(10) 



B. Ordered Polarized Crystals 

The fundamental thermodynamic law determining the 
equations of state for ordered polarized crystals is given 
by 



dw = Tds + E • dP - e : dcr, 



(12) 



wherein w is the enthalpy per unit volume, s is the en- 
tropy per unit volume, T is the temperature, e is the 
crystal strain and cr is the crystal stress. The electric 
adiabatic susceptibility is defined by 



9P 

aE 



(13) 



The tensor describing adiabatic electrostriction co- 
efficients may be defined as^i 



/^^Jfc -\ Qp^ 



(14) 



The piezoelectric tensor is closely related to the elec- 
trostriction tensor via 



( dci 



— l3ijm\mk 



(15) 



When the system is out of thermal equilibrium, the driv- 
ing electric field is 



Erf = E - 



dw 



'dt 



(16) 



wherein r is the relaxation time tensor transport coef- 
ficients which determine the relaxation of the ordered 
polarized system into a state of thermal equilibrium. 



III. RESONANCE DYNAMICS 

Here we shall show how the magnetic intensity 
drives the magnetic resonance equations of motion in 
magnetically ordered systems. Similarly, we shall show 
how the electric field Erf drives the polarization resonance 
equations of motion for polarized ordered systems. 
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A. Gilbert-Landau-Lifshitz Equations 

The driving magnetic intensity determines the torque 
on the magnetic moments according to 



— ^ 7M X Hd 
ot 



(17) 



Employing Eqs.(fTT|) and (|17p . one finds the equations for 
magnetic resonance in the Gilbert form 



dm 

~dt 



7M X 



, (18) 



wherein the Gilbert dimensionless damping tensor a is 
defined as 



a = {jM)t. 



(19) 



One may directly solve the Gilbert equations for the driv- 
ing magnetic intensity according to 



H, + a.(NxH,)=H-(^^ 



(20) 



Egs. pTl) and (l20l) express the magnetic resonance motion 
in the Landau-Lifshitz form. 



B. Landau-Khalatnikov-Tani Equations 

The driving electric field gives rise to a polarization 
response according to 



47r 



(21) 



wherein ujp is the plasma frequency. A simple derivation 
of Eg. (I?!]) may be formulated as follows. In a large vol- 
ume V, the polarization due to charges {zjc} is given 

by 



V 



(22) 



If the driving electric field accelerates the charges accord- 
ing to 



mjVj = ZjeEid, 
then Eq. (j2ip holds true with the plasma frequency 



(23) 



w,^ = 47re^ lim 



V 



= (24) 



wherein Ua is the density of charged particles of type a. 

The polarization resonance equation of motion follows 
from Eqs.ldH) and asii 



'4^\ 5^ 



aP dw{V,s,a) 



dt 



dV 



E. 



(25) 



The electric field E induces the polarization P at reso- 
nant frequencies which are eigenvalues of the tensor 
for which 



? — 1 



(26) 



The decay rates for the polarization oscillations are eigen- 
values of the tensor T for which 



47r 



(27) 



If the decay rates are large on the scale of the the resonant 
frequencies, then the equation of motion is over damped 
so that 



min T j 3> max f2i implies 

j i 

dV , dw{V,s,rT) 
9f + dP 



(28) 



Eg. ([28]) represents the Landau-Khalatnikov equation for 
polarized systems. 



IV. HEATING RATE PER UNIT VOLUME 

Let us here consider the heating rate implicit in relax- 
ation processes. Independently of the details of the mi- 
croscopic mechanism for generating such heat, the rates 
of energy dissipation are entirely determined by t. Ex- 
plicitly, the heating rates per unit volume for magnetiza- 
tion and polarization are given, respectively, by 



dm dm 

qm = —T- ■ T 



at 



dt 



and 



dP dP 

qp = — - T ■ 



dt 



dt 



(29) 



(30) 



Finally, the notion of crystal viscosity r]ijki is introduced 
into elasticity theory^^ via the heating rate per unit vol- 
ume from rates of change in the strain de/dt; It is 



de 



ij deu 



(31) 



Crystal viscosity is employed to describe, among other 
things, sound wave attenuation. Our purpose is to de- 
scribe how heating rates in Eqs. ip^ and ([50)) can be re- 
lated to the heating rate in Ea.f (|3T|) ). This allows us to 
express the transport coefficients t in terms of the crystal 
viscosity. 



A. Relaxation via Magnetostriction 



From the magnetostriction Eq. ([T0| . it follows that 
magnetic relaxation gives rise to a strain 



de. 



dt 



—AijkiNk 
M 



dMi 
dt ' 



(32) 
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and thereby to the heating rate 

4 dM., 



q = J^-g^{AmnqiNq)rimnrsiArsk3Nk)—7^ , (33) 

in virtue of Eq.dST]). Employing Eqs.(l29l) and dSS]), we 
find that the magnetic relaxation transport coefficient in 
the magnetostriction model 



Af2 



(ArskjNk). (34) 



The Gilbert damping tensor follows from Eqs. ([T5|) and 
(l3l as 



close analogy with one another. For the magnetized case, 
the relaxation is driven by the magnetic intensity 
yielding the Gilbert equation of motion^. For the polar- 
ized case, the relaxation is driven by the electric field 
yielding the Tani equation of motion^^. In both cases, 
the relaxation time tensor t is determined by the crystal 
viscosity as derived in the Appendix lA) i.e. in Eqs. (|A3P 
and (jA6[) . The viscosity can be measured independently 
from the magnetic or electrical relaxation by employing 
sound absorption techniques^. 

Appendix A: Kubo formulae 
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Oiij = 'j^i-Ajnnqi'^q^Vranrsi-Arskj-^k)- {^^) 

The central relaxation tensor Eq. ([55l) describes the mag- 
netic relaxation in terms of the magnetostriction coeffi- 
cients and the crystal viscosity. 



From the thermodynamic Eq.ljlOp. the fiuctuations in 
the magnetic intensity are given by magnetostriction, i.e. 



2AijkiNi 
M 



(Al) 



Relaxation via Electrostriction 



Prom the electrostriction Eq. (fT4| , it follows that a time 
varying polarization gives rise to a time varying strain 



de 



dPk 



dt dt 



and thereby to the heating rate. 



dPi 



(36) 



(37) 



Eas. ((XT|) . ^ and © imply 



M2 



{KnnqiNq)Trnnrs{Y,v' ,t){KrskjNk). (A2) 



Employing Eqs. (IA2l) . ^ and ([5]), one finds the central 
result for the magnetic relaxation time tensor; It is 



ma 4 Ctij 

Tij — ji^2 {■^mnqiAfq)Tlfyinrs{Arskj Afk) — ^ji^j ' (-^*^) 



in virtue of Eq.(13T|). Employing Eqs.(|30l) and we 
find that the electric relaxation transport coefficient in 
the electroostriction model 



From the thermodynamic Eq. (|T4|) . the fluctuations in 
the electric intensity are given by electrostriction, i.e. 



(38) 



The central relaxation tensor Eg. ([38)) describes the polar- 
ization relaxation time tensor coefficients in terms of the 
electrostriction coefficients and the crystal viscosity. The 
implications of the electrostriction model for the Landau- 
Khalatnikov equation is to the authors knowledge a new 
result. 



V. CONCLUSIONS 

For ordered polarized and magnetized systems, we 
have developed phenomenological equations of motion in 



AEkir,t) = -/3,j-feAcr.j(r,t). 
Eqs.jXU, ([D) and (P imply 



(A4) 



gff{r, r', t) = PkiiJ^kimn{r, v' ,t)l3„mj- (A5) 



Employing Eqs. dASI) . ^ and ([5]), one finds the central 
result for the electric relaxation time tensor; It is 



^if — PkUVklmnP'. 



mnj • 



(A6) 
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